The Absorption Spectrum of High-Density Stellar Ejecta in the Line-of-Sight 

to Eta Carinae 



in 



T. R. Gull 1 , G. Vieira 1 ' 2 , F. Bruhweiler 1 ' 3 , K.E. Nielsen 1 ' 3 , E. Verner 1 ' 3 ^. Danks 4 

Theodore . R . GullSnasa . gov , GVieiraOst is . gsf c . nasa . gov , FredBOiacs . gsf c . nasa . gov , 
NielsenSstis . gsf c .nasa. gov, KVerner@st is .gsf c .nasa.gov, DanksOngc . com 

o 

^ ■ ABSTR ACT 

> 

Using the high dispersion NUV mode of the Space Telescope Imaging Spectrograph 
(STIS) aboard the Hubble Space Telescope (HST) to observe Eta Carinae, we have re- 
solved and identified over 500 sharp, circumstellar absorption lines of iron-group singly- 
ionized and neutral elements with ~20 velocity components ranging from —146 kms -1 to 
—585 kms -1 . These lines are from transitions originating from ground and metastable 
t^j- \ levels as high as 40,000 cm -1 above ground. The absorbing material is located either in 

dense inhomogeneities in the stellar wind, the warm circumstellar gas immediately in 
the vicinity of Eta Carinae, or within the cooler foreground lobe of the Homunculus. We 
have used classical curve-of-growth analysis to derive atomic level populations for Fe II 
at —146 kms -1 and for Ti II at —513 kms -1 . These populations, plus photoionization 
and statistical equilibrium modeling, provide electron temperatures, T e , densities, uh, 
and constraints on distances from the stellar source, r. For the —146 km s -1 component, 
we derive T e = 6400 K, n H > 10 7 - 10 8 cm -3 , and d w 1300 AU. For the -513 km s -1 
component, we find a much cooler temperature, T e = 760 K, with njr > 10 7 cm -3 , we 
estimate d~ 10,000 AU. The large distances for these two components place the ab- 
sorptions in the vicinity of identifiable ejecta from historical events, not near or in the 
j_j ■ dense wind of 77 Car. Further analysis, in parallel with obtaining improved experimental 
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and theoretical atomic data, is underway to determine what physical mechanisms and 
elemental abundances can explain the large number of strong circumstellar absorption 
features in the spectrum of rj Car. 

Subject headings: stars: individual(Eta Carinae), winds, circumstellar matter - ISM: 
abundances 
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1. Introduction 

Eta Carinae (77 Car) has been an enigma since the Great Outburst in the 1840's (Davidson 
k Humphreys 1999). The 100+ M star, located at a distance ~2300 pc (Meaburn 1999), with 
~10 6 ' 7 L Q and a very large mass loss rate of ~lO _3 M y _1 is surrounded by a dusty, expanding 
double-lobed nebular complex known as the Homunculus. Smith et al. (2003) using infrared 
images estimated that the Homunculus may contain 12 Mq of material, most having been ejected 
since the Great Outburst. Because of its recent activity, proximity, and brightness, r/ Car offers 
a unique opportunity to study mass loss processes of a massive star in considerable detail. Since 
massive Population III stars are thought to have been abundant in the Early Universe (Heger & 
Woosley 2003), the nebulosities of rj Car may provide fundamental clues to how extensive mass 
loss in the last stages of evolution in extremely massive stars, both now and in the Early Universe, 
enriches the interstellar medium in processed elements. 

The observations described herein clearly indicate a rich circumstellar absorption spectrum for 
the ejecta of 77 Car. These narrow, interstellar-like features are seen in roughly 20 components with 
heliocentric velocities from —146 to —584 km s _1 . The spectral lines of this absorption complex 
cannot be interstellar and arise in gas at much higher densities (n ps 10 7 cm~ 3 ) than those of the 
multiple interstellar components seen from +100 to —50 kms -1 toward other stars in the Carina 
association (Walborn et al. 2002). The interstellar lines, typifying densities n < 10 4 cm -3 , are 
formed in the foreground, much lower density H II region gas and interstellar clouds along the 
line-of-sight to n Car. 

2. Observations and Data Reduction 

Observations were obtained with the HST/STIS using the E230H grating and the 0'.'2 x 0'/09 
aperture, providing a spectral resolving power of 114,000 in the near UV. The HST high spatial 
resolution, about 0'.'03, provides critical separation of the stellar emission from dust-scattered 
starlight and nebular emission within the Homunculus. The chaotic nature of the expanding ejecta 
within the Homunculus, the extended wind of the central source and the many emission nebulosities, 
lead to a velocity-smeared spectrum for spatial scales exceeding O'.'l. These two-dimensional STIS 
spectral images were acquired under STIS GTO program 9242 on October 1, 2001 (E230H, 2385 
to 2940 A) and GO program 9083 on January 20, 2002 (E230H, 2885 to 3160 A). Data reduction 
was performed with the CALSTIS software (Lindler 1999) developed at NASA/GSFC. 

3. Analysis 

The STIS spectrum of 77 Car displays sharp circumstellar and interstellar (IS) lines formed 
along the line-of-sight (Figure 1). Walborn et al. (2002) identified multiple IS components in the 
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direction of four association stars (including HDE 303308 and CPD — 59°2603), ranging in velocity 
from +127 to —388 kms" 1 . The strongest IS lines, all arising from levels within the ground 
configuration of the observed species with velocities ranging from +100 to —50 kms -1 , are also 
present in the rj Car spectrum. Unique to the 7] Car spectrum are absorption lines at —146 kms -1 
and the complex circumstellar absorptions between —385 to —585 kms -1 coming from a wide range 
of metastable levels, up to 40,000 cm -1 , in neutral and singly ionized species. Because of the high 
level of excitation, these absorptions imply densities far above what is expected in the normal ISM. 
In Figure 2, absorption line profiles are plotted from —100 to —570 kms -1 for selected lines of 
various species. The slowly varying stellar flux is normalized to facilitate the comparison. Two 
velocity components are readily apparent: one at —146 kms -1 and another at —513 kms -1 . In 
what follows, we will devote our analysis to Fe II and Ti II, respectively, in these two velocity 
components since they are fairly isolated with minimal contamination by absorption from adjacent 
components. 

The two absorption components at —146 and —513 kms -1 are easily distinguishable by their 
very different characteristic widths and by the range in originating energy levels. The relatively 
broad component at —146 kms -1 , observed in lines of Mg II, Cr II, Fe II, Mn II, Ni II and Co II, 
originates from a wide range of energy levels. In Figure 2, the —146 kms -1 component is present 
in Fe II AA2715, 2756, both lines arising from a level 7955 cm -1 above ground. It is significantly 
broader than the multiple, blended components between —385 and —585 kms -1 . This component, 
with identified lines originating from energy levels up to 40,000 cm -1 (nearly 5 eV above ground!), 
is more excited than the higher velocity components. 

The velocity component at —513 kms -1 is observed in lines of both neutral and singly ionized 
species as in Fe I A2995, V II A3111 and Ti II A3088 (Figure 2). While we have identified weak Fe I 
lines originating from levels as high as 7400 cm" 1 , nearly all lines originate from much lower levels. 
The V II absorption lines originate from levels ranging up to 3200 cm -1 (~0.4 eV above ground). 
So far, only upper limits have been found for V II in the interstellar medium (Cardelli 1994; Snow, 
Weiler & Oegerle 1984; Welty et al. 1992), where these results imply vanadium is at least thirty 
times under abundant compared to solar values. In the —513 kms -1 absorption component in the 
spectrum of r) Car, we have unambiguously detected over twenty strong absorption lines of V II. The 
difference in ionization and the relative level populations between the —146 and the —513 kms -1 
velocity components indicate these components originate from regions of very different ionization 
and temperature, and at very different distances from r] Car. 

Extreme care was taken in extracting the information necessary to derive ionic column densi- 
ties in the —146 kms -1 and —513 kms -1 components. Measuring reliable equivalent widths was 
complicated by a) absorption from other components at velocities near the component of inter- 
est, and b) an often highly undulating local continuum due to the complex spectrum produced 
by the underlying photospheric and wind features of the central star. Where the underlying local 
continuum is flat, the equivalent widths and velocity centroids are reliable. But, when the local 
continuum varies sharply with wavelength, the line profile information can be less well defined. 
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Many lines that are obviously present are hopelessly blended or have ill-defined continua; we there- 
for rejected measures of these lines. Fortunately, there are enough strong, well-isolated lines with 
high signal-to- noise to derive useful physical parameters. 

The derived equivalent widths have been analyzed using standard curve-of-growth techniques 
to determine velocity dispersions (or b- values), and column densities for Fe II in the —146 kms -1 
system and Ti II in the —513 kms -1 . The stellar spectrum of rj Car is very complex due to multiple 
P-Cygni-shaped profiles of iron and iron-peak singly-ionized elements. Equivalent widths were 
obtained by a best estimate of the stellar 'continuum'. Estimated errors were obtained assuming 
reasonable limits on the range of possible 'continua'. Thus the error range is much larger than 
what would be expected for photon statistics. For the —146 kms -1 component, Fe II was chosen 
for its large number of spectral lines with equivalent widths lying on the linear portion of the 
curve-of-growth (Table 1 and Figure 3). The —513 kms" 1 component is often blended with other 
components having velocities between —385 and —585 kms -1 , making it more difficult to find 
unblended lines that apply to curve-of-growth analysis (see Figure 2). However, many Ti II, V II 
and Fe I lines are unblended and have proven useful for analysis of the —513 kms -1 component. 
For this discussion we present curves-of-growth for Ti II for the —513 kms -1 component (Table 2 
and Figure 4). 

At this point, we must limit our analysis to those species that have adequate gf-values and 
that have lines which can provide reliable measurements to perform curve-of-growth analysis. We 
are collaborating with several groups to obtain both experimental and theoretical gf-values for the 
lines observed for other species observed in the circumstellar spectrum of rj Car. Consequently, a 
comprehensive, in depth analysis of these features and the species they represent must wait until 
those data are available. 

Variability of the physical conditions is of interest and concern for this ejecta. As noted above, 
the spectra were obtained in two visits separated by three months. Differences in the equivalent 
widths were noticed for Fe II lines (—146 kms -1 component) in the overlap region in the 2900 A 
range. For that reason, we chose to exclude measures of Fe II lines from the January spectrum and 
to use only measured lines up to 2940 A as recorded in the October spectrum. We did not detect 
variations in equivalent widths for the —513 kms -1 absorptions between October 2001 and January 
2002. Hence we used equivalent width measures for Ti II lines from both the October 2001 and 
January 2002 spectra. Variability of the —513 kms" 1 component with time is occurring, but on a 
much longer timescale than for the —146 kms -1 component. The —513 kms -1 component is found 
to be cooler, it is likely further from the star and more shielded by intervening ejecta. This will be 
discussed in a later paper describing changes from the broad maximum beginning in October 2001 
across the 2003.5 Minimum to the early stages of recovery in March 2004. 

We measured equivalent widths for well-isolated Fe II and Ti II spectral lines and used them to 
construct standard curves-of-growth. Column densities were obtained for each energy level using 
transitions with detectable absorption. The curves-of-growth for the three lowest terms of Fe II 
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(Figure 3) are well represented by a curve-of-growth with a b- value of 5.5 kms' 1 , while the Ti II 
curves-of-growth for the two lowest terms (Figure 4) yield b = 2.1 kms -1 . The good agreement 
further attests to the quality of the equivalent width measurements and of the adopted gf- values. 
This is most dramatically seen in the curve-of-growth for the a 4 D term in Fe II and in the a 4 F term 
of Ti II. The very narrow range in the measured velocity centroids (—512.1 ± 1.0 and —145.8 ± 1.0 
kms -1 ), with no systematic variations for lines with larger gf- values and equivalent widths, imply 
negligible contribution from nearby velocities. 

Table 3 presents the derived —146 kms -1 Fe II column densities, plus an estimate for interme- 
diate (14,000 to 21,000 cm" 1 ) energy levels not accessible in this spectral region. The —513 kms 
Ti II has measurable populations for the two lowest configurations (Table 4). The deduced falloff in 
level populations with increasing excitation strongly suggest that contributions from higher energy 
levels to the total column densities are negligible. In this spectral range, we find little evidence for 
Ti I or Ti III for the —513 kms -1 component nor for Fe I or Fe III for the —146 kms' 1 component. 
Thus, it is reasonable to assume that there are only trace amounts of Ti and Fe in other ionization 
states. This allows us to convert derived column densities directly into total elemental column 
densities. 

Because the lowest-lying levels are connected by only parity-forbidden transitions, their popu- 
lations should be well within the conditions approaching LTE for anticipated densities (see Verner 
et al. 1999, 2002). The critical densities for the departure coefficients departing from unity is N cr i t 
= Aji/Cji « 10 6 - 10 8 cm" 3 ( Viotti 1976; see also Figure 3 in Verner et al. 1999). Thus, we 
can use the relative populations of the lowest-lying levels and the Boltzmann Equation to estimate 
the characteristic temperature in the line-forming regions of Fe II and Ti II for the —146 kms -1 
and the —513 kms -1 components, respectively. For the higher energy states (E&), even though 
they are coupled by forbidden transitions, we would expect that their LTE departure coefficients 
to differ from unity. The actual departure coefficients depend strongly on the radiation field and 
cross-sections corresponding to permitted transitions connecting the lower metastable levels with 
higher non- metastable levels. The strong upward radiative rates and subsequent decays can push 
the departure coefficients for the higher energy metastable levels above unity. In actuality, the lower 
level populations (<26,000 cm' 1 ) obey the Boltzman Equation such that the relative populations 
are sensitive only to the electron temperature, T e . Once T e is fixed by the relative populations of 
the lower levels, the departure coefficients of the higher energy levels become a function of density 
and radiation field (Verner et al. 1999). 

We obtained estimates of the electron temperatures, T e , by fitting the Boltzmann equation 
to the populations of the lowest energy levels of Fe II and Ti II. Figures 5 and 6 demonstrate our 
fitting procedure using the Boltzman Equation, Ni/Nq = (gi/ go)e^ Xl ^ x "^ kTe to obtain the electron 
temperature in each component. In deriving these temperatures, we have only used the deduced 
relative level populations (i.e. column density for each energy level) for the lowest energy levels to 
determine temperature. Here, xo denotes the mean energy for levels in the ground terms (See Tables 
3 and 4). For the Fe II in the —146 kms -1 component, our models as presented below indicate 



-6- 



that populations of energy levels are in LTE up through 26,000 cm . As one increases in energy 
beyond 26,000 cm -1 , the models show an increase above unity for the departure coefficients. At 
energies approaching ~40,000 cm" 1 , the departure coefficients are less than five. Above energies of 
«40,000 cm" 1 , the metastable levels disappear and the spontaneous decay rates rapidly depopulate 
the levels. The departure coefficients then approach zero. The net result of our modeling is that 
it appears that our relative populations of Fe II is a function of temperature only. The density 
constraint that we can impose is that the total hydrogen particle density must be log(nn) > 7 to 8. 

For the —146 kms" 1 component, T e = 6400 K, (v t h C rm = 1-5 kms" 1 , compared to a b-value 
of 5.5 kms" 1 ), far below the implied temperature assuming the velocity broadening - deduced 
from the b-value - is purely thermal in origin. This indicates significant turbulence, or possibly 
subcomponents, within each velocity system. In the case of the —513 kms -1 component, the level 
populations of Ti II, which we have assumed the level populations are represented by the Boltzman 
distribution, are consistent with T e = 760 K (vtherm =0.5 kms" 1 , compared to a b-value of 2.1 
kms" 1 ). 

3.1. Modeling the Velocity Components 

We have applied photoionization simulations, incorporating statistical equilibrium calculations 
for a 371-level model for Fe + to produce a grid of models representing the absorbing gas of the 
the —146 and the —513 kms" 1 components. Similar simulations have been used successfully to 
model the Fe II and [Fe II] emission arising in the Weigelt Blobs B & D (Verner et al. 2002) 
- ejecta of rj Car associated with the 1890's event (Smith et al 2004). The models in our grid 
are calculated based upon measured column densities, temperatures, and relative level populations 
derived from the curve-of-growth analysis. From our previous work (Verner et al. 2002), we adopt a 
stellar luminosity of lO^ergs" 1 with T e g- = 15,000 K. We approximate the stellar flux distribution 
by a Kurucz (1979, 1988) model atmosphere. We assume that the ejecta have undergone CNO- 
processing, where carbon and oxygen are at 0.01 solar abundance, with nitrogen at lOx solar 
abundance (Verner, Bruhweiler, & Gull 2004; also see Hillier et al. 2001). Given our poor overall 
knowledge about dust formation processes in the Eta Carinae environment, uncertainties in gas- 
phase fraction, dust particle size and its chemical compositions, the effects of dust are not included 
into our current calculations. However, the geometry of the Little Homunculus indicates that a hot, 
wind-blown cavity is between the star and the Little Homunculus (Ishibashi et al. 2003, Nielsen 
et al. in prep). We thus expect little or no dust between the central source and the —146 kms" 1 
component. 

The temperature of the low velocity (—146 kms" 1 ) component is very close to that previously 
derived for the Weigelt Blobs B and D (Verner et al. 2002). Thus, it is reasonable to expect 
an ionization for this component to be similar to the Weigelt Blobs. Unfortunately, our spectral 
data do not provide a means to estimate the (Fe/H) abundance ratio. On the other hand our 
emission line study (Verner, Bruhweiler and Gull 2004) as well as the investigation by Hillier et al. 
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(2001) of the stellar spectrum show that (Fe/H) ranges from 0.2x (Fe/H) so i ar in the BD Blobs to 
2x(Fe/H) so i ar , respectively. Therefore given the uncertainties, we have adopted the (Fe/H) to be 
(Fe/H) so i ar in our calculations. 

The most accurate measurements of the relative Fe II level populations (or column densities) 
are for the second (a 4 F) and third (a 4 D) configurations. We have constructed the theoretical level 
population ratio and scaled the other Fe II level populations derived from the curves-of-growth to 
this ratio. The Fe II relative energy level population ratio (R=n(a 4 D)/n(a 4 F)) is plotted in Figure 
7 for a range of hydrogen densities from 10 5 to 10 8 cm -3 and a range of stellar distance between 
10 15 - 5 and 10 16 - 5 cm. The best models suggest the temperature range is from 5700 K to 7300 K, well 
bracketing the Fe II-inferred temperature of 6400 K (Figure 5). The predicted Fe II populations 
for the a 4 F and a 4 D terms are found to be in LTE. In Figure 8 the relative Fe II population levels 
are compared to the predicted level populations for njj = 10 7 and 10 8 cm -3 at distances 10 16 and 
10 16 ' 5 cm. From Figure 7, the normalized (unity) ratio indicates the stellar distance for the —146 
kms" 1 to be 10 16 3 cm, or «1300 AU. 

Independently, we repeated calculations for the —513 kms -1 component. We measured only 
column densities for Ti II from the two lowest configurations as indicated above. Consequently, we 
cannot derive total hydrogen densities for this absorbing region. However, given the low deduced 
T e and low ionization of this component, one must conclude that if it were at a distance comparable 
or less than the —146 kms -1 component, it would necessarily have a higher density. Otherwise, 
the ionization and temperature would be much higher than the T e = 760 K found. From a purely 
qualitative argument, we conclude the —513 kms -1 component, unless it has exceptionally high 
densities, is at a distance much farther than 1300 AU. 

For convenience, we adopt a nn >10 7 cm" 3 , similar to that found for the 146 kms -1 component. 
Again, we assume the same stellar luminosity, but use the temperature derived from the Ti II energy 
populations, T c =760 K. Our photoionization modeling then places the —513 kms -1 component 
eight times more distant, or r«10,000 AU from the stellar source. 

Based upon our modeling, we present the deduced physical conditions for the —146 kms^ 1 and 
the —513 kms -1 components in Table 5. Of course, these results are somewhat dependent upon 
the gaseous elemental depletions for Fe (—146 kms -1 ) and the Ti (—513 kms -1 ). We have also 
ignored possible shielding in these calculations. For example, if other components, lying between 
the central stellar source and the —146 kms -1 and —513 kms -1 , produce significant attenuation 
of the T) Car radiation field, then this would reduce the implied distance for that component. 

4. Interpretation 

The geometry of the Homunculus, based upon STIS CCD long aperture spectra of nebular 
emission and absorption lines originating from the interior surface of the expanding debris, has 
been determined to be a bipolar lobed structure at an inclination angle of ~41° degrees with 
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respect to the plane of the sky (Davidson et al. 2001). Spatial variations along the slit are resolved 
at the O'.'l level, a projected scale of about 230 AU at the distance of rj Car. The line-of-sight, 
centered on rj Car, passes almost tangentially through the foreground lobe of the Homunculus, and 
likely intersects clumps of corresponding size. The Homunculus is a relatively thin shell. However, 
the analysis of Davidson et al (2001) used the emission of [Fe II] and [Ni II] lines to describe what 
appeared to be a relatively diffuse, but thin interior surface of the Homunculus. N. Smith (private 
communication) comments that the infrared spectroscopy of the Homunculus is consistent with a 
shell about one-tenth the thickness of the distance from r] Car, and with a relatively diffuse interior 
skin and very thin outer skin. 

The —513 kms -1 component is likely part of the Homunculus in the line-of-sight. Based upon 
its relatively small b-value, it may be associated with the thin outer skin. Likewise, the detected 
components with velocities between —385 and —585 kms -1 are evidence for clumps within the 
relatively thick wall of the Homunculus. Interior to the Homunculus is a large cavity, filled with 
a hot, very low density gas, too faint to be detected by current observations. Still interior is the 
Little Homunculus (Ishibashi et al. 2003) that is a thin, ionized shell with measured velocities in 
the vicinity of —146 kms -1 . Expansion velocities and proper motions of both the Homunculus and 
the Little Homunculus are consistent with ejection during the major outburst in the 1840's and 
the lesser outburst in the 1890's. The central region, immediately surrounding rj Car, is thought 
to be a wind-blown cavity, with an extended skirt whose plane is perpendicular to the axis of the 
two lobes. In this picture, the two major velocity components seen in the STIS NUV spectra arise 
in the gas ejected during each of the two outbursts. The other velocity components, intermediate 
in velocity space, likely are components of the outer Homunculus, distributed in that one-tenth 
thickness, but projected along a long line of sight due to geometry. We note in passing that these 
multiple components are photoionized (or excited) by rj Car, and that during the minimum, very 
significant changes occur. Data have been obtained during the 2003.5 minimum with the STIS, 
analysis on the temporal variability is ongoing. 

The presence of dust and its effects on the excitation of these ejecta components appears to be 
surprisingly little. However, given the hot, low-density cavities bounded by the Little Homunculus 
and by the interior wall of the Homunculus, little dust intervenes between the star and the Little 
Homunculus, and in turn between the Little Homunculus and the outer Homunculus. The bulk of 
the dust is likely contained in the clumps of ejecta in line of sight. Hence we might expect to see 
some dust contribution in the thin wall of the Little Homunculus (—146 kms -1 component) and 
considerably more dust contributions in the thicker wall of the Homunculus. Indeed dust is present 
in the Homunculus and is the scatterer for the stellar radiation seen reflected of the bipolar lobes. 
Smith et al (2003) demonstrated that the Balmer alpha (Ha) P-Cygni line profile changed across 
the 5.54-year cycle consistent with the terminal wind velocity increasing to about 1200 kms -1 
during the maximum and dropping to about 600 kms -1 during the minimum. These changes in 
wind velocity are interpreted to originate from the high latitudes of the central source and that 
the dust in the relatively flat polar regions of the Homunculus has an optical depth of order unity, 
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therefore causing a single scattering event between the star and the observer. Given the very 
lumpy appearance of the flat polar region and the relatively linear striations of the walls of the 
Homunculus, dust appears to be clumpy in distribution in the outer structure of the Homunculus. 
The simplest interpretation is that most of the dust in line of sight resides in the outer wall of the 
Homunculus and that some may lie in the interior thinner wall of the Little Homunculus. Certainly 
the presence of Fe II absorption lines in the Little Homunculus and Fe I absorption lines in the 
Homunculus indicates that much of the FUV radiation is stopped by the Little Homunculus. Some 
FUV radiation down to the Lyman limit does get through the wall as the STIS FUV and FUSE 
spectra show complex, but relatively continuous UV radiation to well below 1000 A originating from 
a nearly point-like source. However no locally emitted Balmer line emission is detected from the 
Homunculus, but near-red [Fe II] and [Ni II] emission lines are present and were used by Davidson 
et al. (2001) to trace the geometrical structure of the Homunculus. The emission line velocities 
of the Little Homunculus and the Homunculus in the line of sight are in good agreement with the 
absorption velocities we see. Thus, the simplest interpretation of these absorption components are 
that they are associated with the expanding ejecta, not clumps within the much hotter wind within 
a few hundred AU of the star. 

Many of the elements identified in these velocity components are refractory, with condensation 
temperatures below 1300 K. In a later paper we plan to address relative abundances, but initial 
inspection indicates that elemental depletion is low and not following a Spitzer-Routly law (Routly 
& Spitzer 1952). In particular the presence of V II and Ti II at —513 kms -1 is remarkable. We 
speculate that the combination of low C and O abundances, the expanding ejecta and/or the 
modulated wind may inhibit normal grain formation in the r\ Car system. 

5. Conclusions 

The UV STIS spectrum of 77 Car reveals complex multiple, narrow absorption features for 
many neutral and singly-ionized species, spanning a velocity range from —146 kms -1 to —585 
kms -1 . Two narrow velocity components (—146 and —513 kms -1 ) appear to originate from shells 
or clumps in the foreground lobe of the Homunculus. These clumps are heated by the very bright 
UV flux of 77 Car. The presence of ionic transitions up to energies of 40,000 cm -1 is, in part, 
due to the high density in ejecta giving rise to the observed very strong circumstellar absorptions. 
The photoionization modeling for the —146 and the —513 km s -1 components suggest that the 
corresponding absorptions arise at large distances from the central stellar source. We conclude that 
the absorption is formed in dense gas located at different positions within the bipolar structure of 
the Homunculus. The —146 kms -1 component is formed at a distance of ~1300 AU, which leads 
to a spectrum consisting mainly of lines from singly-ionized iron-group ions. The —146 kms -1 
component is likely formed in the wall of the Little Homunculus (Ishibashi et al. 2003). Although 
the conclusion is less firm, the cooler —513 kms -1 component, seen in the lines of neutral and singly 
ionized elements, appears to be formed at a much larger distance from the star (~10,000AU). This 
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component is likely formed in the wall of the Homunculus. Further analysis of the STIS dataset for 
rj Car, including changes that occur during the 2003.5 minimum will be presented in subsequent 
papers. 

This paper is based upon observations made with the NASA/ESA Hubble Space telescope, 
obtained at the Space Telescope Science Institute, which is operated by the Association of Univer- 
sities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. Funding was provided 
through the STIS Guaranteed Time Observations (GTO). The research of EV has been supported 
through STIS GTO funding and through NSF - 0206150 to CUA. 
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Table 1. Fe II lines used in the —146 kms 1 component analysis 



Lower Level 


Aiab a 


Vel. b 


Jlow 


Eiow d 


w A c 


log(9/) f 


Ref. 




(A) 


(kms- 1 ) 




(cm- 1 ) 


(A) 






( J D)4s a 6 D 


















ZOoD.DOU 


1 A Q Q 


A K 

4.0 


U.UUU 




(1 in 

— o.iy 


r 1V1 W 




2600.173 


— 145.9 


4.5 


0.000 


0.264 


0.35 


FMW 




2599.147 


— 145.2 


3.5 


384.790 


0.165 


—0.10 


TVTV 111 T 




2399.973 


— 144.8 


2.5 


667.683 


0.184 


—0.15 


FMW 




2407.394 


— 144.9 


1.5 


862.613 


0.176 


—0.25 


TV TV 111 T 

(MW 




2614.605 


— 142.9 


1.5 


862.613 


0.180 


—0.39 


TV TV 111 T 

(MW 




2621.191 


— 146.0 


1.5 


862.613 


0.050 


— 1.83 


FMW 




2414.045 


— 145.2 


0.5 


977.053 


0.165 


—0.43 


FMW 




2622.452 


— 146.1 


0.5 


977.053 


0.119 


— 1.00 


TVTV 111 T 

(MW 


d art 


















2451.842 


— 146.4 


4.5 


1872.567 


0.031 


—2.07 


K88 




2484.951 


— 146.5 


4.5 


1872.567 


0.047 


—1.98 


K88 




2392.206 


— 146.7 


3.5 


2430.097 


0.056 


—1.64 


FMW 




2505.969 


— 147.0 


3.5 


2430.097 


0.022 


—2.22 


K88 




2512.125 


-146.6 


3.5 


2430.097 


0.018 


-2.53 


K88 


( 5 D)4s a 4 D 


















2563.304 


-144.7 


3.5 


7955.299 


0.156 


-0.05 


FMW 




2693.633 


-146.5 


3.5 


7955.299 


0.026 


-2.09 


FMW 




2715.218 


-146.2 


3.5 


7955.299 


0.256 


-0.44 


FMW 




2740.358 


-145.8 


3.5 


7955.299 


0.203 


0.24 


FMW 




2756.551 


-144.8 


3.5 


7955.299 


0.224 


0.38 


FMW 




2881.601 


-145.5 


3.5 


7955.299 


0.041 


-1.66 


FMW 




2927.441 


-146.9 


3.5 


7955.299 


0.085 


-1.23 


FMW 




2564.245 


-146.0 


2.5 


8391.938 


0.151 


-0.29 


FMW 




2869.716 


-147.0 


2.5 


8391.938 


0.012 


-2.29 


FMW 




2940.367 


-146.7 


2.5 


8391.938 


0.008 


-2.75 


FMW 




2567.683 


-145.3 


1.5 


8680.454 


0.098 


-0.65 


FMW 




2731.542 


-145.7 


1.5 


8680.454 


0.082 


-0.95 


FMW 




2769.753 


-146.4 


1.5 


8680.454 


0.047 


-1.51 


FMW 




2762.629 


-147.1 


0.5 


8846.768 


0.068 


-1.28 


FMW 



a Laboratory wavelengths 

b Measured velocity of absorption line 

c J-values for the lower energy state 

d Energy of the lower state 

c Measured equivalent width 

f log(gf) values from references 

References. — (FMW) - Fuhr J. R., Martin G. A, and Wiese W. L. 1988, J. Phys. 
Chem. Ref. Data 17, Suppl. 4, (K88) - Kurucz R. L. 1988, In: McNally M. (ed.) Trans. 
IAU, XXB, Kluvcr, Dordrecht, p. 168 
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Table 2. Ti II Lines used in —513 kms 1 Component Analysis 



Lower Level 




Vel. 


T c 

" low 




Wa 


l °s (a j) 






(Km s ) 




(cm ) 


(A) 




( 3 F)4s a 4 F 
















3058.282 


-513.1 


1.5 


0.00 


0.016 


-1.78 




3067.238 


-513.5 


1.5 


0.00 


0.053 


-0.71 




3073.863 


-512.5 


1.5 


0.00 


0.070 


-0.32 




3122.515 


-514.4 


1.5 


0.00 


0.006 


-2.36 




3148.948 


-513.0 


1.5 


0.00 


0.043 


-1.22 




3060.634 


-513.9 


2.5 


94.10 


0.012 


-1.57 




3067.109 


-513.3 


2.5 


94.10 


0.052 


-0.59 




3076.117 


-512.5 


2.5 


94.10 


0.080 


-0.12 




3131.706 


-513.0 


2.5 


94.10 


0.026 


-1.19 




3158.307 


-513.7 


2.5 


94.10 


0.008 


-2.17 




3073.000 


-512.9 


3.5 


225.73 


0.046 


-0.62 




3079.539 


-512.6 


3.5 


225.73 


0.081 


0.06 




3144.665 


-512.4 


3.5 


225.73 


0.019 


-1.30 




3088.922 


-512.1 


4.5 


393.44 


0.085 


0.23 


d 3 b 4 F 
















2536.632 


-512.1 


1.5 


908.02 


0.009 


-1.03 




3155.106 


-512.2 


1.5 


908.02 


0.016 


-1.15 




3162.117 


-512.1 


1.5 


908.02 


0.023 


-0.69 




2535.381 


-512.9 


2.5 


983.89 


0.012 


-0.94 




3153.163 


-510.7 


2.5 


983.89 


0.012 


-1.06 




3162.684 


-512.0 


2.5 


983.89 


0.032 


-0.55 




2532.013 


-512.7 


3.5 


1087.32 


0.012 


-0.67 




3156.582 


-510.9 


3.5 


1087.32 


0.008 


-1.17 




3163.481 


-512.1 


3.5 


1087.32 


0.033 


-0.38 




2526.362 


-512.8 


4.5 


1215.84 


0.014 


-0.51 



a Laboratory wavelengths S. Johansson (private communication) 
b Measured velocity of absorption line 
c J-values for the lower energy state 
d Energy of the lower state 
e Measured equivalent width 

f log(gf) values from Pickering, Thorne & Perez (2001,2002) 
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Table 3. Column densities from the —146 kms Fe II curve-of-growth analysis. 



Conf. 


Emean (™ 1 ) 1 


ni (10 13 cm- 2 ) 




416 


1 fiQ Q 


3d 7 a 4 F 


2417 


200.0 


4sa 4 D 


8321 


54.56 


a 4 P 2 


18007 


59 3 


4sb 4 P 


21421 


0.83 


4sa 4 H 


21459 


3.85 


4sb 4 F 


22808 


2.34 


4sa 4 G 


25761 


0.95 


4sb 2 P 


26169 


0.51 


4sb 2 H 


26253 


0.99 


4s a 2 F 


27445 


0.86 


4sb 2 G 


30556 


0.42 


4sb 4 D 


31421 


0.68 


4sa 2 I 


32892 




4sc 2 G 


33482 


0.18 


4sc 2 D 


38184 


0.28 


N^n, 4 




495 



1 Emean— X)i ( E igi)/ Si 9i 

2 Missing energy states intermediate in energy between 
4sa 4 D and 3d 7 b 4 P 

3 Estimated column density based on interpolation 

4 For completeness, we also note that column densities for 
energy states a 6 S, b 2 F, b 2 D, and a 2 S were also not mea- 
sured but as at high energy levels would contribute a negligible 
amount to the total population. 



Table 4. Column densities from the —513 kms Ti II curve-of-growth analysis. 





Conf. 


Emean (cm 1 ) 1 


n, (10 13 cm~ 2 ) 




4s a 4 F 


225 


16 




3d 3 b 4 F 


1085 


3 


N = J2i n i 






19 



1 Emean — Y^i (Eigi)/ J2i 9i 



-15- 



Table 5. Physical parameters for the —146 and —513 km s 1 components 





— 146 kms 1 


-513 kms" 1 


b- value 


5.5 kms - 1 


2.1 kms" 1 


Derived T e : 


6400 K 


760 K 


Ionization 


Singly Ionized Species 


Neutral & 






Singly Ionized Species 


Hydrogen density, n e : 


>10 7 cm~ 3 


>10 7 cm- 3 


Distance from star 


~1300 AU 


~10,000 AU 



MO 



-11 



8-10 



4-10 



2-10 



Eta Car-A 
HDE 303308 
CPD -59 2603 




-600 




II 1 \ ,' II 

II 



:i 




-400 



-200 



Velocity (km s"') 



Fig. 1— Spectral region of Mg I A2853 for rj Car, HDE 303308 and CPD -59° 2603, all members 
of the Carinae association. Strong intervening IS absorption velocities are +100 to —50 kms -1 
(Walborn et al. 2002). Unique to the r] Car spectrum are the prominent absorption line at —146 
kms -1 and the multiple, narrow absorption lines between —385 and —585 kms -1 . The rj Car 
spectrum includes a slowly varying P-Cygni profile in addition to a series of narrower absorption 
components between —146 and —585 kms -1 . 
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Fig. 2. — Comparison of circumstellar absorption line profiles for selected observed transitions. 
These normalized spectral profiles are for transitions from a range of lower levels for several neutral 
and singly-ionized elements. The —513 kms -1 component is the narrowest, yet absorption lines 
originate from levels up to 7400 cm _1 above ground. Note that many Fe I, V II and Ti II lines are 
unique to the —513 kms -1 component. The —146 kms -1 component is from singly- ionized, iron- 
peak and other elements, includes transitions from even higher energy levels and is characteristically 
broader than the lines of the —513 kms -1 component. The complicated velocity structure between 
—380 to —513 kms -1 is noticeable for the two Fe II lines. The two Fe II lines also have a very 
broad wind component, starting from —250 kms^ 1 and extending into the portion of the spectrum 
of the narrow line components. Spatial resolution of these structures demonstrate that the wind 
feature is very distinct from the narrow velocity components. 
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a 6 D 



a 4 F 



a 4 D 



b=10.0 
b=H.() 




b = 5.5 km/s 

N = 1.7 x 10 15 cm" 2 
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b=8.0 
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b = 5.5 km/s 

N = 2.0 x 10 15 cm" 2 
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b=6.0 
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N = 5.5 x 10 14 cm" 2 



3 

log (gtt) 



Fig. 3. — The —146 kms -1 Component curves-of-growth plots for the lowest three terms of Fe II. 
Measured lines were selected to minimize blending with other velocity components, to be measurable 
against the relatively complex stellar spectrum. The relatively large uncertainty bars are due to 
ability to determine the continuum level of rj Car, and not photon statistics. Note that all three 
curves-of-growth are consistent with a b- value of 5.5 kms -1 . 
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b=4.0 




log (gfX) 



Fig. 4. — The —513 kms Component curves-of-growth plots for the lowest two terms of Ti II. The 
— 513 kms -1 component is much narrower in width and is seen in lines originating from significantly 
lower energy levels that seen either in the P-Cygni wind of 77 Car or in the multiple absorption 
components similar in excitation to the —146 kms -1 component. No lines originating from higher 
energy levels were found. Thus the populations of these two terms are found to be the major 
contributors to the total amount of Ti II. 
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Fig. 5. — Fe II population for each energy level of the —146 kms component. The linear slope 
corresponds to 6370 K. Note that the population for each energy level is normalized by the statistical 
weight. 
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Fig. 6. — Ti II population for each energy level of the —513 kms component. The linear slope 
corresponds to 760 K. 
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Fig. 7. — The —146 kms component: hydrogen density versus distance plot. The normalized ratio 
of the Fe II column densities for the a 4 F and a 4 D terms, R=(a 4 D/a 4 F) t h e ory/(a 4 D/a 4 F) b S ervation, 
was used to determine the optimal density/distance value. A value of unity for R indicates a match 
to the derived column density. At a distance of 10 161 to 10 16 ' 4 cm, the hydrogen column density 
would be between 10 to 10 8 cm~ 3 . 
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Fig. 8. — The observationally-derived and predicted Fe II relative level populations for the —146 
kms -1 component. The best fit lies between the two density values, possibly closer to nn=10 8 
cm -3 . The theoretical and observed relative level populations roughly approximate a straight line 
for the three cases shown. The slope reflects the electron temperature, T e , for each dataset. Most 
of the scatter is due to differing statistical weights for each level. 



